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Theoretically, some works have proposed the hadronic resonances Ec(2800) and Ac(2940)^ to be 
S-wave DN and D*N molecular candidates, respectively. In the framework of QCD sum rules, 
we investigate that whether Ec(2800) and Ac(2940)"^ could be explained as the S-wave DA'' state 
with = i and the S-wave D* N state with — | , respectively. Technically, contributions 
of operators up to dimension 12 are included in the operator product expansion (OPE). The final 
resuhs are 3.52 ± 0.36 GeV and 3.59 ± 0.41 GeV for the S-wave DN state of = i" and the S- 
I wave D* N state of = | , respectively. They are somewhat bigger than the experimental data of 

Ec(2800) and Ac(2940)^, respectively. These results do not support the statements that Ec(2800) 
. could be the S-wave DN molecule with = 5 and Ac(2940)+ be the S -wave D* N molecule 



(N 



with J = f , respectively. As byproducts, masses for their bottom partners are predicted to be 



^ ! 6.97 ± 0.55 GeV for the S-wave BN state of = i and 6.98 ± 0.56 GeV for the S-wave B*N 

state of 
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PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 



^ i I. INTRODUCTION 

Ph. 

Oh' In the past several years, many new excited charmed baryonic states have been discovered experimentally. 
(— I , For example, Belle Collaboration observed an isotriplet of new states Sc(2800) decaying into A+tt, and they 
tentatively identified the quantum numbers of these states as = | fl. In particular, the neutral state 



Sc(2800)° was possibly confirmed in the Sc(2800)''p channel by Babar Collaboration [3]. However, 

K*" ' the measured mass 2846 ± 8 ± 10 MeV for Ec(2800)° is 3cr higher (assuming Gaussian statistics) than the 
, Belle's measured value, and Babar indicated that there is weak evidence that the excited S° observed by 
' them is = i. Moreover, Babar collaboration reported the observation of a new charmed state Ac(2940)+ 
^ ' decaying to D°p with a mass of 2939.8 ± 1.3 ± 1.0 MeV and an intrinsic width of 17.5 ± 5.2 ± 5.9 MeV 

. Subsequently, Belle Collaboration confirmed it in the Ac(2940)+ — > Ec(2455)°'"'""'"7r+'~ decay and measured 
^ ; its mass and width to be 2938.0 ± 1.3^1'^ MeV and 13^5!^ MeV, respectively Q. 

y—i ' The experimental observations have triggered theorists' great interest in understanding the inner struc- 
^ ■ tures of I]c(2800) and Ac(2940)+. One direct way of theoretical studies grounds on the assignments of 
, them as conventional charmed baryons. From a relativized potential model prediction, masses of I]c(2800) 
*^ ' and Ac (2940)+ are close to theoretical values of S* with = |~ or |~ and A* with = |" or I"*", 
_Cd_' respectively [5|- In the relativistic quark-diquark picture, Ebert et al. suggested Sc(2800) as one of the 
orbital (IP) excitations of the Ec with = ^ , § ' o'" I ' ^^'i proposed Ac(2940)+ as the first radial 
excitation of Ec with = |^ 6]. Strong decays of Ec(280G) and Ac(2940)"'" as charmed baryons were 
analyzed in heavy bary on chiral perturbation theory [tI, with the ^Pq model Q, and using the chiral quark 
model [91]. In Ref. [l0|, Garcilazo et al. indicated that Ec(2800) would correspond to an orbital excitation 
with = i or I and Ac (2940)"*" may constitute the second orbital excitation of the Ac by the Faddeev 
method. In a mass loaded flux tube model, Chen et al. suggested that Ac(2940)+ could be the orbitally 
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excited A+ with = | ll|. He et al. evaluated the production rate of Ac(294G)+ as a charmed baryon 



Another different way bases on the assumption that Ec (2800) and Ac (2940)+ are some molecular can- 



didates. Lutz et al. interpreted Ec(2800) as a chiral molecule [13|. In Ref. [14], Ec(2800) was deciphered 
as a dynamically generated resonance with a dominant DN configuration. Dong et al. pursued a possible 
hadronic molecule interpretation of Ec(2800) as a bound state of the charmed D meson and the nucleon 
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N, since I]c(2800) isotriplet states are very close to respective DN thresholds 15|. They chose several 
possible quantum number assignments of Sc(2800) as = ^ and | , and concluded that Sc ~^ A^tt 
decay widths are consistent with current data for the — and — | assignments. Here — \ 
corresponds to an S'-wave DN configuration, = and = |^ represent a P-wave, and — | 
has a relative D-wave in the DN system. Coming down to Ac (2940)+, it was firstly proposed to be a 
S'-wave D*^p molecular state with = ^ in Ref. [16|, because its mass is just a few MeV below the 
D*^p threshold. From an effective Lagrangian approach, the strong two-body decay of Ac(2940)"'" was 
studied under the = ^ and D*N molecular assignments, and it showed that — i should be 
ruled out Later, the radiative and strong three-body decays of Ac(2940)+ were also researched in the 



D*N molecule picture with = [18|,llS!- He et al. systematically studied the interaction between 



D* and N, and concluded that the D*N systems may behave as = and |^ states [20|. In Ref. 



2l|, Recio et al. found a possible molecular candidate for the Ac (2940)"'" in the | channel. Ortega et al. 
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studied Ac(2940)"'" as a D*N molecule with J = ^ in a constituent quark model, and gained a mass 
which agrees with the experimental data [22]. 

Although various interpretations to Ec(2800) and Ac(2940)''" were put forward, at present their underly- 
ing structures are still unclear, which means that it is interesting and significative to make more theoretical 
efforts to reveal their properties. Therefore, we devote to studying that whether Ec(2800) and Ac(2940)"'" 
could be the S-wave DN state with = i and the S'-wave D*N state with = | , respectively. 
QCD is widely believed nowadays to be a true theory of strong interactions. At high energy, the effective 
coupling constant of the quark-gluon interaction becomes small because of asymptotic freedom and the 
interaction can be treated perturbatively. On the other hand, quark interaction within hadrons is strong 
since it binds quarks into unseparable pairs. Thus, low energy QCD involves a regime where it is futile to 
attempt perturbative calculations, and the strong interaction dynamics of hadronic systems is governed by 
nonperturbative QCD effects completely. There are still many questions on nonperturbative QCD remain 
unanswered or realized only at a qualitative level since one's absence of knowledge on QCD confinement 
effects. Therefore, it is quite difficult to calculate the hadron spectrum from QCD first principles. The 
method of QCD sum rules |23j . developed by Shifman, Vainshtein, and Zakharov, represents an attempt 
to bridge the gap between the perturbative and nonperturbative sectors by employing the language of 
dispersion relations. It is well known for the advantages of this method: instead of a model-dependent 
treatment in terms of constituent quarks, hadrons are represented by their interpolating quark currents 
and the interactions of quark-gluon currents with QCD vacuum fields critically depend on the quantum 
numbers (spin-parity, flavor content) of these currents. The QCD sum rule method is a nonperturbative 
formulation firmly based on the first principle of QCD, which has become a widely used working tool in 
hadron phenomenology. The mere fact that the seminal paper on QCD sum rules have already been cited 
more than 4000 times reflects its vigorousness. It has been successfully applied to conventional mesons 



and baryons (for reviews see 2J-|27| and references therein) and multiquark states (e.g. see [2 




particular, many theoretical practitioners began to study light pentaquark states in Refs. 
heavy pentaquark systems in Ref. [3l| . 

In this work, we construct interpolating currents for the S-wave DN state with = i and the 
S-wave D*N state with = | , and then investigate from QCD sum rules that whether I]c(2800) and 
Ac(2940)+ could be the S-wave = \~ DN state and the S-wave = |~ D*N state, respectively. 
The rest of the paper is organized as follow. We derive QCD sum rules for molecular states in Sec. HIl 
with similar techniques as our previous works on heavy baryons 32| as well as molecular states 33|. The 
numerical analysis and discussions are presented in Sec. IIIIl and masses of D'^*'>N and B(*)iV molecular 
states are extracted out. Sec. IIVI contributes to the conclusions. 
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II. QCD SUM RULES FOR MESON-NUCLEON MOLECULAR STATES 

A. interpolating currents 

The basic point of QCD sum rule calculations is to construct a proper interpolating current to represent 
the studied state. At the present time, currents for S'-wave D'^*^N or B'-*^N molecular states can be 
constructed from corresponding currents of D*-*^ or i?^*) mesons and A'^ nucleons to form meson-baryon 
confi gur ations of fields. In full theory, interpolating currents for Z?(*' and S'*) mesons can be found in 



Ref. [34:], and currents for nucleons have been listed in Ref. 35 1. Therefore, we build the following forms 
of currents: 

J = {Qc'n5q,'){eabcqI''C^^qh5rq'3), (1) 
for ^-wave DN or BN molecular state with spin parity — i , and 

jP = {Q,'7''qc')ieabcqI''Cj^qh5rql), (2) 

for S'-wave D*N or B*N molecular state with spin parity = ^ ■ Here Q is heavy quark c or b, and qi, 
q2, as well as (73 denote light quarks u and/or d. The index T means matrix transposition, C is the charge 
conjugation matrix, with a, 6, c and c' as color indices. 

B. QCD sum rules for meson-nucleon molecular states 

The QCD sum rules for DN and BN molecular states with = i are constructed from the two-point 
correlator 



n(q2) = j j d4^e^?-(0|T[j-(^)j(0)]|0). (3) 
Lorcntz covariance hints that the correlator has the form 

n(g2)^^ni(g2)+n2(q2). (4) 
Phenomenologically, the correlator can be expressed as 



^2 i-Mu r ^ JlmnP^^'^(g) -H mnP^^^(,s) 



where is the mass of the hadronic resonance, sq is the threshold parameter, and Xh gives the coupling 
of the current to the hadron {0\j\H) — Xhu{p, s). In the OPE side, the correlator can be written as 

IJrnl S~q J S-q'' 

After equating the two sides for IV(q^), assuming quark-hadron duality, making a Borel transform, and 
transferring the continuum contribution to the OPE side, the sum rules can be written as 
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where indicates the Borel parameter. 

There is some difference on deriving of mass sum rules for D*N and B*N molecular states with — | . 
One can start from the two-point correlator 

W^iq^) = I J d^xe^'^-^{0\T[f{x)r{0)m- (9) 

Lorentz covariance implies that the two-point correlator in Eq. ([9]) has the form 

U'^q') = -gP^m,{q')+Il,{q')] + (10) 

where the ellipse denotes other Lorentz structures which acquire contributions from both J = ^ and J = | . 
The tensor structures g'''^4 ^^'^ 5''^ contributed only by the J = ^ hadrons. The phenomenological 
side of UP'^{q'^) can be expressed as 

where Xh gives the coupling the coupling of the hadronic state to the current as {0\j'^\H) = Xnv^lq, s). 
Here, v''{q, s) is the Rarita-Schwinger spinor for = I .In the OPE side, one can write the correlator 



as 



U''^q') = -9 



( \ 



"d,^+n^ond(^2)l ^12^ 



Equating the two sides for n'''^(g^), assuming quark-hadron duality, and making a Borel transform, the 
sum rules can be expressed as 

y2^^-MllM- ^ / dspi(s)e-^/^^'+i3n$=0'^d, (13) 

To eliminate the hadron coupling constant \h in sum rules jZ]), dH), ([T3|) . and (fT4|) . one can take the 
derivatives of sum rules with respect to 1/M^, divide the results by themselves to get 



Ml 



Ml 



d.pi(s)se-«/^Vd/d(-^)Bn?°nd(^)|/|^^^'' d5pi(.)e-^/*^VBn?ond(^)| (^5) 
d.p2(.)se-^/^^VdM-^)Bn^ond(,)|/|^"' dspM^-^l''^ 



C. spectral densities 

The spectral density is given by the correlator's imaginary part 

p,(s) = ilmnpPE(s), J ^1,2. (17) 

In the concrete OPE calculation, one works at leading order in as and considers condensates up to di- 
mension 12. To keep the heavy-quark mass finite, one can use the momentum-space expression for the 
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heavy-quark propagator 34 1 
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{<7r,\{i> + mo) + (|5 + mg)crKA] 



(18) 



+ Toff / *-^7<5*-^<5e'-'en 
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The hght-quark part of the correlator can be calculated in the coordinate space employing the light-quark 
propagator 



Sab{x) 



iSab ^ _ mqSab 



{gqa ■ Gq) 



ix Sab I - /-y \ / ^ab /_ \ / 2/~i2\ 



jmq{gqa ■ Gq):/; 



(19) 



3.26 ^^■^ 27.32 -^^^-^ ^'^ 210-33 

which is then Fourier-transformed to the momentum space in D dimension. The resulting light-quark part 
is combined with the heavy-quark part before it is dimensionally regularized at Z? = 4. 
The spectral densities can be listed as 

/ X pert. N , (qq) , N , (qq)^ / \ , (gqcr-Gq) , x , (g^G^), n , (g^G'^}, x , (qq)" f \ , iqq) (gq<y-Gq) , x 



_piS<q^-Gq)(gq.-Gq)^^^^^^iMq)(a^G^)^^^^^pm^9^^^^^ ,^^^2. 



(20) 



Because many terms of p2{s) are proportional to light quarks' masses and approximate to zero, we merely 
present spectral densities resulted from Iii{q'^). Concretely, they read 
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for the S'-wave DN or BN state with = i , and 
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for the S'-wave D*N or i3*A^ state with = | . The lower hmit of integration is given by A = itiq/s. 



III. NUMERICAL ANALYSIS AND DISCUSSIONS 

We perform numerical analysis of the sum rule ([15]) to extract mass values of studied states. One could 
take input parameters as = 1.23 ± 0.05 GeV, = 4.24 ± 0.06 GeV, (qq) = -(0.23 ± 0.03)^ GeV^, 
(gqa-Gq) = (qq), = 0.8±0.1 GeV^ (g^G^) = 0.88 GeV^ and (g^G^) = 0.045 GeV^ 0. To choose 
proper work windows for the threshold Sq and the Borel parameter M'^, one could consider two rules in 
the standard QCD approach: on one hand, the perturbative contribution should be larger than condensate 
contributions to have a good convergence in the OPE side; on the other hand, the pole contribution should 
be larger than the continuum state contributions to satisfy "single-pole+continuum states" hypothesis in 
the phenomenological side. Besides the above two restrictions, the threshold parameter ^io should not be 
taken arbitrarily since it characterizes the beginning of continuum states. 

However, there may be some difficulty to find a conventional work window critically sat isfy ing both of 



the above two rules here, which have been discussed in detail for some other cases, e.g. 36l-l39l|. The main 
reason is that some condensate contributions are very large, making the standard OPE convergence (i.e. 
perturbative contribution at least larger than each condensate contribution) to happen only at very large 
values of M'^. To obtain some useful hadronic information from QCD sum rules, one could try releasing the 
critical OPE convergence criterion and consider the ratio of perturbative contribution to the "total OPE 
contribution" (the sum of perturbative and other condensate contributions calculated) instead of the ratio 
of perturbative contribution to each condensate contribution. Actually, the present case is not too bad, 
there are denumerable important condensate contributions and some of them could cancel out each other 
to some extent since they have different signs. Most of other high dimension condensate contributions 
are very small. All these factors bring that the OPE convergence is still under control at relatively low 
values of M'^. It is also very important that the ratio of perturbative contribution to the "total OPE 
contribution" has been examined and does not change much even including some condensate contributions 
higher than dimension 12. 

The dependence on Borel parameter A'P for masses of S'-wave DN and BN states with = ^ are 
shown in FIG. 1, for which continuum thresholds are taken as ^/sq = 3.2 ~ 3.4 GeV and ^/sq = 6.7 ~ 
6.9 GeV, respectively. From the Borel curves, one can visually see there indeed exist stable plateaus and one 
could expect the OPE convergence is under control. Thus, we choose some transition range 3.0 ^ 4.0 GeV^ 
as a compromise Borel window and arrive at 3.52 ± 0.16 GeV for DN state. Considering the uncertainty 
rooting in the variation of quark masses and condensates, we gain 3.52 ± 0.16 ± 0.20 GeV (the first error 
reflects the uncertainty due to variation of y/so and M'^, and the second error resulted from the variation of 
QCD parameters) or concisely 3.52 ±0.36 GeV for the S'-wave DN state with = i , which is somewhat 
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higher than the experimental value of Ec(2800) even considering the result's uncertainty. This does not 
support the statement that £^(2800) could be a S-wave DN molecule with = ^ , which agrees with 
the conclusion that A^tt decay widths are not compatible with data for the = ^ molecular 

assignment of I]c(2800) in Ref. 15]. For the 5- wave BN state with = i , we choose some transition 
range 4.0 ^ 5.0 GeV^ as a compromise Borel window and arrive at 6.97 ± 0.34 GeV. Considering the 
uncertainty rooting in the variation of quark masses and condensates, we gain 6.97 ± 0.34 ± 0.21 GeV or 
6.97 ± 0.55 GeV for the S'-wave BN state with = i". 
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FIG. 1: Masses of S'-wave DN and BN states with — ^ as a function of Af ^ from sum rule (|15|l are shown. 
The continuum thresholds are taken as ^/so — 3.2 ~ 3.4 GeV and ^/so = 6.7 ~ 6.9 GeV, respectively. 
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FIG. 2: Masses of 5- wave D* N and B* N states with 
The continuum thresholds are taken as ^/so — 3.3 



as a function of M from sum rule l|15p are shown. 
3.5 GeV and = 6.7 ~ 6.9 GeV, respectively. 



Masses of 5- wave D*N and B*N states with = | as a function of A/^ from sum rule (IT5|) are 
shown in FIG. 2. Graphically, one can see there have stable plateaus for Borel curves. Similarly, we 
choose some transition range 2.5 ^ 3.5 GeV^ as a compromise Borel window for D*N state, and arrive at 
3.59 ±0.22 GeV. Considering the uncertainty rooting in the variation of quark masses and condensates, we 
obtain 3.59 ± 0.22 ± 0.19GeV or 3.59 ± 0.41 GeV in a concise form for the S'-wave D*N state with = 
I , which is bigger than the experimental data of Ac(2940)+ even taking into account the uncertainty. 
Therefore, it does not support the statement that Ac(2940)+ could be the S'-wave D*N molecule with 



J' 



We choose a compromise Borel window 4.0 ~ 5.0 GeV^ and take V^o = 6.7 ~ 6.9 GeV for 
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the S-wa.ve B*N state with = |~. In the work windows, we obtain 6.98 ± 0.35 GeV for B*N state. 
Varying input values of quark masses and condensates, we attain 6.98 ±0.35 ±0.21 GeV or 6.98 ±0.56 GeV 
for the S'-wave B*N state with J-^ = |". 

IV. CONCLUSIONS 

In some theoretical approaches, the hadronic resonances Ec(2800) and Ac(2940)+ have been suggested 
to be S'-wave DN and D*N molecular states, respectively. From QCD sum rules, we investigate that 
whether 1:^(2800) and Ac(2940)+ could be the S'-wave DN state with = i~ and the D*N state with 
— I , respectively. In the OPE calculation, contributions of operators up to dimension 12 are included 
to examine that its convergence is still under control. The final result for the 5-wave DN state of — ^ 
is 3.52 ± 0.36 GeV, which is somewhat bigger than the experimental value of Sc(2800) even considering 
the result's uncertainty. The numerical result for the S'-wave D*N state of = | is 3.59 ± 0.41 GeV, 
which is higher than the experimental data of Ac(2940)"''. These results do not support the statements 
that 1:^(2800) and Ac(2940)+ could be the S'-wave DN molecule with = and the D*N molecule 
with = I , respectively. By the way, we have also studied corresponding bottom counterparts and 
predicted their masses to be 6.97 ± 0.55 GeV for the S-wave BN state with = i and 6.98 ± 0.56 GeV 
for the B*N state with = | , which could be searched in further experiments. 
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